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Abstract. Equilibrium systems which exhibit a phase transition can be studied by 
investigating the complex zeros of the partition function. This method, pioneered 
by Yang and Lee, has been widely used in equilibrium statistical physics. We show 
that an analogous treatment is possible for a nonequilibrium phase transition into 
an absorbing state. By investigating the complex zeros of the survival probability 
of directed percolation processes we demonstrate that the zeros provide information 
about universal properties. Moreover we identify certain non-trivial points where the 
survival probability for bond percolation can be computed exactly. 



PACS numbers: 02.50.-r,64.60.Ak,05.50.+q 
1. Introduction 

The investigation of nonequilibrium systems is of great importance since most 
phenomena in nature take place under nonequilibrium conditions. In this field, 
as in equilibrium physics, phase transitions are particularly interesting. However, 
when dealing with nonequilibrium systems one cannot utilize such a well-established 
theoretical framework as in equilibrium statistical mechanics. Therefore it is interesting 
to investigate which concepts of equilibrium physics can be transfered to nonequilibrium 
systems. 

An important concept of equilibrium statistical physics is the theory of Yang and 
Lee [1] for the emergence of nonanalytic behaviour at phase transtitions. Only recently 
these ideas have been applied to an integrable nonequilibrium model [2]. In the present 
work we use similar techniques to investigate directed percolation (DP) [3] as a paradigm 
for nonintegrable systems far from equilibrium. We show that it is possible to transfer 
the ideas of Yang and Lee to DP and that this method does indeed provide information 
about universal properties of the phase transition. 

% To whom correspondence should be addressed (daimer@comphys.uni-duisburg.de) 
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1.1. Yang and Lee's theory in equilibrium statistical physics 

We shall shortly sketch the ideas of complex zeros of the partition function in equilibrium 
statistical mechanics. Taking the Ising model without an external magnetic field as an 
example, we consider a system of N spins (s, = ±1, i — 1, . . . , N) on a lattice where 
the spins interact with their nearest neighbours, denoted by < i,j >. The total energy 
7i({sj}) of a configuration {sj} of the spins is given by 7~t({si}) = —J^2 <i j > SjSj. The 
system is said to be in thermal equilibrium at temperature T when the probability 
distribution P({sj}) to find the system in configuration {sj} is given by the stationary 
Gibbs ensemble P({sj}) oc e _/3W ^ Si ^, (3 = kvT, where A; B is Boltzmann's constant. In 
that case the statistical properties of the system are fully determined by the canonical 
partition function 

Z N (T) = J2^' m{s ' }) = J2 xr({St}) with (!) 

{Si} {Si} 

x = ef" > , r({ Sl }) = s i s j ■ ( 2 ) 

<i,j> 

For finite system size Zn can be expressed as the ratio of two polynomials in x with 
integer coefficients. The free energy per lattice site is obtained from Z N through 

The hallmark of phase transitions is the appearance of nonanalyticities. Since in 
Equation (3) singularities occur as zeros of the partition function Z^, these are good 
candidates for indicating a phase transition. However, Z^ has no zeros for a positive 
real temperature so that there is no phase transition for finite N. This led Yang and 
Lee to study the complex zeros of Z N ^, which in this case correspond to complex 
temperatures, and examine their behaviour as the system size grows. They argued that 
in the thermodynamic limit iV — > oo in (3), if the model exhibits a phase transition, 
an ever increasing number of complex zeros accumulate in the vicinity of a point on 
the positive real axis, which corresponds to a physical temperature. The zeros thereby 
induce singularities at the critical point in the thermodynamic limit, thus explaining the 
crossover to nonanalyticity at the transition in the limit of an infinite system size. Such 
a behaviour has been observed for many equilibrium models. Moreover, several features 
of the distribution of zeros have been related to universal properties of the system under 
consideration [4]. 

1.2. Directed percolation, a process away from equilibrium 

Since nonequilibrium systems do not obey the stationary Gibbs ensemble they cannot 
be handled using the partition function of equilibrium statistical physics. Instead they 
are described by the time-dependent probability distribution P t ({sj}) which has to 

For a field-driven phase transition these zeros are denoted as Yang-Lee zeros while in the case of a 
temperature-driven phase transition the term Fisher zeros is used. 
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be derived from the master equation. In most cases this cannot be done exactly. 
Nevertheless, the concept of universality, well-known from equilibrium statistical 
physics, proves to be suitable for nonequilibrium phase transitions as well. In this 
context models exhibiting a transition from a fluctuating active phase to a nonffuctuating 
inactive phase (absorbing state) have been extensively studied [5-8]. These models 
are used to describe spreading processes as e.g. forest fires, where a spreading agent 
can either spread over the entire system or die out after some time. The two phases 
are seperated by a nonequilibrium phase transition. The most prominent universality 
class of transitions into an absorbing state is that of directed percolation. Directed 
percolation is an anisotropic variant of ordinary percolation in which activity can only 
percolate along a given direction in space. Regarding this direction as a temporal 
degree of freedom, DP can be interpreted as a dynamical process. Directed percolation 
emerges in a variety of physical problems ranging from catalytic reactions on surfaces [9] 
to spatio-temporal intermittency in magnetic fluids [10]. 

1.3. Realizations of DP 

Simple realizations of DP are directed bond (DPb) and directed site percolation (DPs) 
on a tilted square lattice (see Figure 1). In directed bond percolation the bonds are 
conducting with probability p and non-conducting with probability 1 — p. In this model 
sites at time t > are activated by directed paths of conducting bonds, originating 
from active sites at time t — 0. A cluster consists of all sites that are connected by such 
paths of conducting bonds to active sites in the initial state. In directed site percolation 
on the other hand all bonds are conducting while the sites themselves can be either 
permeable (p) or blocked (1—p). Activity can spread from permeable site to permeable 
site. A cluster is formed by permeable sites that are connected to active sites at time 
t = by a directed path of bonds that only connect permeable sites. 

The order parameter which characterizes the phase transition is the probability 
P(oo) that a randomly chosen site belongs to an infinite cluster. For p > p c this 
probability is finite whereas it vanishes for p < p c . Close to the phase transition P(oo) 
is known to vanish algebraically as P(oo) ~ (p — Pc) 13 , where (3 is an universal critical 
exponent. In addition, the DP process is characterized by a spatial correlation length 
(perpendicular to time) and a temporal correlation length £||. As p approaches p c 
these length scales are known to diverge as 



where x is the position and z = v\\/v±_ is the so-called dynamic exponent. Numerical 
estimates for the values of p c and the critical exponents in 1 + 1 dimensions are 
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"time" t 




"space" coordinate j 

Figure 1. Lattice geometry of DP on a tilted square lattice in 1+1 dimensions. In 
both cases (DPb and DPs) sites can be eiter active (sj tt = 1) or inactive (sj tt = 0). 
In DPb activity can percolate forward in time through bonds (solid lines) which are 
conducting. In DPs acitvity can percolate through bonds that connect permeable sites. 

p c = 0.6447001(1), = 0.27649(4), is ± = 1.096844(14) and v\\ = 1.733825(15) [11] 
(bond percolation) and p c = 0.705489(4) [12] (site percolation). 

Although DP can be defined and easily simulated, it is one of the very few systems 
for which - even in one spatial dimension - no analytical solution is known, suggesting 
that DP is a non-integrable process. In fact, the values of the percolation threshold 
and the critical exponents are not simple numerical fractions but seem to be irrational 
instead. 

1.4. Yang and Lee zeros for directed percolation 

To apply the idea of Yang-Lee zeros to DP, we consider the order parameter in a finite 
system as a function of the percolation probability p in the complex plane. This can be 
done by studying the time-dependent survival probability P{t), which is defined as the 
probability that a cluster generated in a single site at time t = survives up to time 
t (or even longer). For a finite system P(t) can be expressed as a polynomial in p (see 
below). Note that lim^oo P(t) and the order parameter P(oo) coincide. 

The survival probability and the partition function (1) show a similar behaviour 
in many respects. For finite systems P{t) and Z N do not have relevant zeros in the 
physical region of the control parameter < p < 1 and < T < 00 although the phase 
transition is marked by a vanishing P(t) and at the critical points in the limit of 
infinite systems. At those points both functions exhibit nonanalytic behaviour. The 
zeros of P(t) and Zn are generated by polynomials with real integer coefficients, thus 
the zeros come in complex-conjugate pairs. 

According to Yang and Lee the increasing system size is accompanied by an 
approach of some of the complex zeros of Z N to the physical region of the control 
parameter, while the accumulation point marks the critical point. We show that the 
complex zeros of P(t) approach the critical value p c on 'trajectories' for increasing time 
t and that the distance between the zeros and the critical point is related to the critical 




bond percolation site percolation 



Figure 2. Example of clusters which survive until time t — 3. Acitve sites are shown 
as black circles, conducting bonds as solid lines. Left: Dashed lines denote the hull, 
irrelevant bonds are not shown. The corresponding weight to the survival probability 
P(3) in (6) is p 4 (l — p) 4 . Right: White patterned circles belong to the hull. Simple 
white circles are irrelevant for the survival of the cluster. The corresponding weight to 
P(3) is p\l -pf. 



exponent for the temporal correlation length. 
1.5. Determining the survival probability 

In directed bond (site) percolation the survival probability P(t) is given by the sum over 
the weights of all possible configurations of bonds (sites) for which the process survives 
at least up to time t. Each conducting bond (permeable site) contributes to the weight 
with a factor p, while each non-conducting bond (blocked site) contributes with a factor 
1 — p. However, the states of those bonds (sites) which do not touch the actual cluster 
are irrelevant as they do not contribute to the survival of the cluster. Therefore, it is 
sufficient to consider the sum over all possible clusters C of bonds (sites) connected to the 
origin. Each cluster is weighted by the contributions of the conducting bonds (permeable 
sites) belonging to the cluster and the non-conducting bonds (blocked sites) belonging 
to its hull. Roughly speaking, the hull surrounds the cluster. More precisely the hull 
of a cluster is the set of non-conducting bonds (blocked sites) that would contribute 
to the cluster if they were conducting (permeable) (see Figure 2). Thus the survival 
probability can be expressed as 

P(f) = 5> B ( 1 -*0 m > (6) 
c 

where the sum runs over all clusters reaching the horizontal row at time t. For each 
cluster n denotes the number of its bonds (sites), while m is the number of bonds (sites) 
belonging to its hull. Note that in this sense the hull does not include bonds connecting 
sites at time t and t + 1 (DPb) or sites at time t + 1 (DPs) since the cluster may survive 
even longer. Summing up all weights in Equation (6) one obtains a polynomial. As can 
be verified, the first few polynomials for the survival probability in the directed bond 
percolation process are 
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P(0) = 1 (7) 
P(l) = 2p-p 2 

P(2) = 4p 2 - 2p 3 - 4p 4 + 4p 5 - p 6 

P(3) = 8p 3 - 4p 4 - 10p 5 - 3p 6 + 18p 7 + 5p 8 - 30p 9 + 24p 10 - 8p u + p 12 
P(4) = 16p 4 - 8p 5 - 24p 6 - 8p 7 + 6p 8 + 84p 9 - 29p 10 - 62p u - 120p 12 

+244j9 13 + 75p 14 - 470p 15 + 495p 16 - 268p 17 + 83p 18 - 14p 19 + p 20 . 

As t increases the number of cluster configurations grows rapidly, leading to complicated 
polynomials with very large coefficients (e.g. for t — 15 the largest coefficient for bond 
percolation is of order 10 44 ). 

The paper is structured as follows. In Section 2 we investigate if the distribution 
of zeros is totally dependent on microscopic details of the underlying model or whether 
it also contains universal features. We ask whether the distribution is related to some 
of the critical exponents which characterize the phase transition of directed percolation. 
Section 3 presents certain non-trivial points where the polynomials P(t) for bond 
percolation can be calculated exactly for all times. Section 4 ends with a conclusion. 
In Appendix A we show that the value of the survival probability can be computed by 
Monte Carlo simulations even if p is a complex number. The first coefficients of the 
polynomials P(t) for bond percolation are discussed in Appendix B. 



2. Universal properties 

The distribution of the zeros of P(t) (from t — 8 to t — 15) in the neighbourhood of the 
critical point is shown in Figure 3. Away from the critical point the appearance of the 
distributions for DPb and DPs is quite different. This implies that in this region the 




3 



Re(p) 



Figure 3. Zeros of the survival probability P(t) in the complex plane. 
Left: Directed bond percolation. Right: Directed site percolation. 
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Figure 4. Zeros approaching the critical point. Three 'trajectories' are enumerated 
and shown as dashed lines. Left: Directed bond percolation. Right: Directed site 
percolation. 



distribution of zeros depends strongly on microscopic details and hence is non-universal. 
However, there is also a general feature which can be observed in both cases. The 
innermost zeros approach the critical point on 'trajectories' as t increases (see Figure 4). 
We applied a standard Bulirsch-Stoer (BST) acceleration algorithm [13] to the set of 
zeros of each enumerated 'trajectory' of Figure 4 to determine lim^ooRe (p zero (t)) and 
lim^oo Im (p zero (t)). The results are listed in Table I. Although we calculated the 
zeros only for small systems (until t — 15) the extrapolants accord fairly well with the 
numerical values of the percolation threshold" 1 " p c = 0.6447001(1) (bond percolation) and 
p c = 0.705489(4) (site percolation). Thus a similar scenario as observed by Yang and 
Lee for equilibrium phase transitions proves to be suitable for the nonequilibrium phase 
transition of DP. As time increases the zeros of the survival probability P(t) approach 
the real axis between p = and p = 1. The accumulation point is the critical point. 

So far we have shown that the zeros of P(t) provide information about the existence 
of the phase transition and about the critical value p c . However, the value of p c is non- 
universal and depends on the particular realization of DP. On the contrary universal 
features are independent of the microscopic details of the underlying process. On each 

+ The convergence for site percolation is slower than for bond percolation since the order of the 
polynomials P(t) for DPs is smaller than for DPb. 



trajectory 


Re (pTZ°J 




trajectory 


Re (pTZ°J 


Im(p£°J 


bond 1 


0.64472(1) 


0.0001(1) 


site 1 


0.70547(3) 


0.0001(3) 


bond 2 


0.6445(2) 


0.008(1) 


site 2 


0.709(4) 


0.006(1) 


bond 3 


0.6470(4) 


0.051(7) 


site 3 


0.712(1) 


0.01(1) 



Table 1. Bulirsch-Stoer extrapolants for the zeros of the 'trajectories' from Figure 4 
for both bond and site percolation in 1+1 dimensions. 
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trajectory 


1/^11 


trajectory 




bond 1 


0.57675(3) 


site 1 


0.5765(7) 


bond 2 


0.575(2) 


site 2 


0.5771(2) 


bond 3 


0.576(4) 


site 3 


0.5731(2) 



Table 2. Bulirsch-Stoer extrapolants for the exponent On each 'trajectory' in 

Figure 4 the distance d(t) — \p zcr °(t) — p c \ is assumed to decrease as d(t) <~ i -1 /"!!. 

'trajectory' of Figure 4 we calculated the distance between the zeros and the critical 
point d(t) = \p zcro (t) — p c \ with the values of p c as given in section 1.3. According to 
Equation (4) and simple scaling arguments we expect d(t) to decrease as d(t) ~ 
Application of the BST algorithm yields results which support this claim, as shown 
in Table 2. The extrapolants are in fairly good agreement with the numerical value 
of 1/V| = 0.57676, even for small systems. This means that universal properties of 
the phase transition of directed percolation are encoded in the complex zeros of the 
survival probability P(t). Figure 5 shows the used data and the asymptotic power law 
d(t) ~ t~ x l v w with the numerical value of v\\. 




time t time t 



Figure 5. The distance from the critical point d = \p zcm — p c \ versus time t. For 
comparison the power law d <~ t~ x l v \\ is shown, with v\\ — 1.7338. Left: Directed bond 
percolation. Right: Directed site percolation. 

3. Exact results 

In this section we will address a particularly surprising observation, namely the existence 
of certain points on the real axis where the polynomials P(t) for bond percolation can 
be solved exactly for all values of t. Beside the trivial points p = (where P(t) = 8 t ,o) 
and p — 1 (where P(t) = 1) we find a t-independent zero at p = 2 and, even more 
surprisingly, a very simple solution if p is equal to one of the Golden Ratios (1 ± v / 5)/2. 
The Golden Ratios are the roots of the quadratic equation p 2 — p + 1 and play an 
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Figure 6. Example for different weights for P(t) and R(t). Left: The contribution of 
the shown cluster to P(3) in (6) is p 4 (l — p) 2 . Right: The contribution of the cluster to 
R(3) in (9) is p 4 (l — p) 6 since it is required that active sites at t — 3 are not connected 
to sites at t — 4 by active bonds. 



important role not only in number theory [15] but also in other fields ranging from 
chaotic systems [16] to arts [17]. Although these special points are located outside 
the physically accessible region < p < 1, their existence may be helpful for further 
investigations of the polynomials P(t). 

3.1. Time-independent zero at p = 2 

For p = 2 and t > 1 all polynomials P(t) vanish identically. This can be shown as 
follows. Let us consider the probability R(t) that a cluster dies out at time t, i.e., the 
row at time t is the last row reached by a cluster. Obviously R(t) is related to the 
survival probability by 



Clearly, R(t) can be expressed as a weighted sum over the same set of clusters as in (6). 
However, in the present case the weights differ from those in Equation (6) by the number 
of non-conducting bonds in the clusters hull between t and t + 1 since it is now required 
that all sites at time t+1 are inactive (see Figure 6). This means that R(t) can be 
expressed as 



where n, m and C have the same meaning as in (6) and k is the number of active sites 
in the horizontal row at time t. Obviously, for p = 2 the additional factor (1 — p) 2k 
drops out so that P(t) = R(t) for all values of t. Moreover, R(0) = P{0) = 1 (for 
p = 2). Combining these results with Equation (8) we arrive at P(t) = for t > 0, 
which completes the proof. 

3.2. Exact solution for p at the Golden Ratio 

For p = (l±y / 5)/2 we find that the survival probability 'oscillates' between two different 
values, namely 



R(t) =P(t)-P(t+l). 



(8) 




(9) 



c 




if t is even 



if t is odd. 



(10) 



Yang-Lee zeros for a nonequilibrium phase transition 



10 



To prove this result, we first verify that (10) is indeed satisfied for t = and t — 1. 
Then we show that 

P{t) = P{t -2) fort>2 and p=(l±V5)/2. (11) 

However, instead of analyzing the survival probability directly, it turns out to be more 
convenient to consider the complementary probability Q(t) = 1 — Pit) that a cluster 
does not survive until time t. Obviously, Q(t) is the sum over the weights of all clusters 
which do not reach the horizontal row at time t, i.e., we impose the boundary condition 
so,t = sij = . . . = St,t = 0. Depending on the states of the two sites so,t-i and so,t-2 
at the left edge of the clusters, this set of clusters may be separated into three different 
subsets, namely, 



(a) a subset where So,t-i — s o,t-2 = 1, 

(b) a subset where So,t-i = and So,t-2 = 1, and 

(c) a subset where So,t-i = So,t-2 = 0. 



Next we show that the weights of the clusters in the subsets (a) and (b) cancel each 
other. To this end we note that the weighted sum Q(t) over all clusters in subsets 
(a) and (ft) may be decomposed into two independent factors Q(t) = Q1Q2, where Qi 
depends only on the state of the three bonds between the sites So,t-2> so,t-i> So,t> and Si it 
(inside the box in Figure 7), while Q2 accounts for all other relevant bonds. Obviously, 
the first factor is given by 

Qf ] =p{i-P)\ Q? = i-p, (12) 

while Q2 takes the same value in both subsets. Thus, if p is given by the Golden Ratio, 
we obtain + Qf^ = Qi(t) = and therefore the weights of subsets (a) and (b) 
cancel each other. Consequently, all remaining contributions to Q(t) come from the 
clusters in subset (c) where the sites So,t-i an d «o,t-2 are inactive. Now we can iterate 
this procedure by successively considering the sites Sj,t-i and from the left to the 



t-2 
t-1 

t O D 



7^ 



000 

(a) 



O 

o o 



000 
(b) 



o 

o 
o o 



000 

(c) 



Figure 7. Decomposition of the configurational sum of non-surviving clusters into 
three subsets. Open and closed bonds are denoted by solid (dashed) lines. Bonds which 
are not shown may be cither open or closed. The box includes all bonds contributing 
to the factor Q\ while all other bonds contribute to the factor Q^. The proof shows 
that the configurations in (a) and (b) cancel each other so that only the configurations 
of (c) contribute to Q. Iterating the procedure by shifting the box to the right, it can 
be shown that all sites at t — 2 and t — 1 have to be zero. 
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right, where j — 1 . . .t — 2. In this way it can be shown that all these sites have to 
be inactive as well. Therefore, the only surviving contributions are those in which the 
entire row of sites at t — 2 is inactive, implying that Q(t) = Q(t — 2). The proof of 
Equation (10) then follows by induction. 

4. Conclusions 

In this paper we have investigated the applicability of Yang-Lee theory to the 
nonequilibrium phase transition of directed percolation. In Section 1 we developed 
an idea of how to transfer the concepts of Yang and Lee to the phase transition of 
DP by studying the complex zeros of the survival probability. In Section 2 we showed 
that a similar scenario as observed for equilibrium phase transitions is also suitable 
for DP, namely complex zeros approaching the critical point for increasing system size. 
Moreover, we could extract the value of the critical exponent for the temporal correlation 
length from the distribution of zeros. Hence the zeros encode univerals properties of 
the phase transition of DP. Section 3 presented exact results on the survival probability 
of bond percolation which may be helpful for further investigations. More precisely, 
we proved that there are certain non-trivial values of p for which the values of the 
polynomials P(t) can be calculated exactly for all times. 
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Appendix A. MC-Simulations 

Here we demonstrate that the value of the survival probability for a complex percolation 
probability p as well as the polynomials P(t) are accessible by computer simulations. 

Let us consider the definition of the survival probability as a sum over configurations 
of surviving clusters in Equation (6). Formally this expression can be rewritten as 



where n and m denote the number of bonds of the cluster and of its hull respectively. 
Therefore, instead of simulating the system at a given p, we may simulate it using a 
different percolation probability q reweighting each cluster by the factor f n ,m- While 
q £ (0, 1) still has to be a real number, p is no longer restricted to be real, it can be any 
complex number. Using this reweighting technique it is in principle possible to access 
the entire complex plane by numerical simulations. For any finite t such a simulation 




(A.l) 
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is stable and will converge to the correct result. But even for small t the convergence 
time can be very long, limiting the range of applications. 

Using the reweighting technique it is also possible to approximate the coefficients 
of the polynomial P(t) = ^ k aj,p k by considering p as a free parameter and expanding 
the term (1 — p) m in Equation (A.l). The coefficients are then given by 



where q G (0, 1) is again a free parameter and 



if < k — n < m 



C{m,k-n) = { (*-«)Km-k-n)i ~"- n (A.3) 
[ otherwise . 

However, in most cases the direct construction of the polynomials using symbolic algebra 
turns out to be more efficient. 



Appendix B. First coefficients of the polynomials P(t) in the limit t — > oo 

As can be seen in Equation (7) the first non-vanishing coefficient of the polynomial 
for the survival probability P(t) = ^2 n a n p n for DPb is always a power of 2. This 
contribution corresponds to the configurational weight of a surviving path without 
branches. The purpose of this appendix is to point out that even the following 
coefficients can be computed exactly, provived that t is large enough. More specifically, 
we conjecture that 

a t+k = 2 t - k - 1 q k (t), (B.l) 



where 



m=0 



is a finite polynomial in t. Using the explicit expressions for P(t) up to t — 15 we find 
that the first eight polynomials read 

<fo(t) = 2 
gi(i) = -2 
q 2 (t) = -4-2* 
q 3 (t) = -2t 

q 4 (t) = - 16 + 3t + t 2 (B.3) 
q 5 (t) = 32 + t + M 2 

q 6 (t)= - 232+yt + 3t 2 -^ 3 

Q7 (t) = 808- 51 + 5t 2 - |t 3 

We conjecture that the leading coefficient of these polynomials is given by 

(fc) J 2(-l)"/M! if fc even 

I 2(fc- 1)(-1) M /M! if fc odd . 1 ' 
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This implies that the first no n- vanishing coefficients of the polynomial P(t) grow in 
such a way that their limit lim^oo for fixed k is well-defined. Summing up these 
contributions we obtain for the survival probability 

P(t) ~ (2pf J2 t&BL (i + p{j - 1/2)) (B.5) 

3=0 J ' 

= h y 2p) t e- p2t/ \A-2p-pH) . 

Physically this expression corresponds to loop-free graphs and serves as an 
approximation close to p — 0. The convergence radius of Equation (B.5) determined 
from 2p = e p2/4 is \p\ < 0.53744. Note that none of the non-trivial roots computed up 
to t — 15 lies inside this radius. We conjecture that this might be true for any t. 
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